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Abstract We investigate pseudogap phenomena originating from pairing fluctu- 
ations in the BCS-BEC crossover regime of a two-dimensional Fermi gas in a 
harmonic trap. Including pairing fluctuations within a T-matrix theory and effects 
of a trap within the local density approximation, we calculate the local density 
of states (LDOS) at the superfluid phase transition temperature T c . In the weak- 
coupling regime, we show that the pseudogap already appears in LDOS around the 
trap center. The spatial region where the pseudogap can be seen in LDOS becomes 
wider for a strong pairing interaction. We also discuss how the pseudogap affects 
the spectrum of the photoemission-type experiment developed by JILA group. 
PACS numbers: 03.75.Hh, 05.30.Fk, 67.85.Lm. 



1 Introduction 



Ultracold Fermi atom gases are now widely recognized as a powerful quantum 
system to study many-body physics in strongly interacting fermion systems^ 2 -. 
One typical example is the BCS (Bardeen-Cooper-Schrieffer)-BEC (Bose-Einstein 
condensation) crossover phenomenon realized in 40 K and 6 Li Fermi gases, where 
the character of superfluidity continuously changes from the BCS-type to the 
BEC-type, with increasing the strength of a pairing interaction associated with 
. a Feshbach resonance 12 - 3 . Using this phenomenon, we can study Fermi superfiu- 

^ ids and Bose superfiuids in a unified manner. The intermediated coupling regime, 

which is sometimes referred to as the crossover region, has attracted much at- 
tention, because various anomalous phenomena originating from strong pairing 
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fluctuations are expected there. Indeed, anomalous single-particle excitation spec- 
tra have been recently observed in a 40 K Fermi gas by using a photoemission-type 
experiment 4 - 5 . Theoretically, it has been pointed ou t 6 i 7 > 8 i 9 that this anomaly can 
be quantitatively understood as a pseudogap phenomenon associated with strong 
pairing fluctuations. 

Although the origin of the observed anomalous single-particle excitation spec- 
tra is still in debate, the above mentioned pseudogap scenario is a strong candi- 
date. In this regard, to assess the validity of this idea, a two-dimensional Fermi gas 
would be useful, because stronger pairing fluctuations than the three-dimensional 
case are expected. Recently, the photoemission-type experiment has been done 
for a quasi-two-dimensional 40 K Fermi gas 10 . In this experiment, in addition to 
the so-called back-bending curve, the split of the spectrum is also observed in the 
crossover region, which is quite different from the three-dimensional result (where 
a single broad peak is only observed). 

In this paper, we investigate single -particle excitations in the BCS-BEC crossover 
regime of a two-dimensional Fermi gas. Although the recent experiment deals 
with a quasi-two dimensional gas 10 , we ignore the z-direction, for simplicity. To 
examine the validity of the pseudogap scenario proposed in the three-dimensional 
case, we include strong-coupling pairing fluctuations within a T-matrix theory in 
the normal state. We also include effects of a harmonic trap within the local den- 
sity approximation (LDA). Using this combined T -matrix theory with LDA, we 
calculate the local density of states (LDOS), as well as the photoemission spec- 
tra, at the superfluid phase transition temperature T c . We note that this framework 
has succeeded in quantitatively explaining the observed photoemission spectra in 
a three-dimensional 40 K Fermi gas 8 ' 9 . We also note that, in the presence of a trap, 
the present T-matrix theory with LDA gives a finite T c , as in the BEC case of 
a trapped Bose gas 12 . In the absence of trap, the superfluid phase transition in a 
uniform two-dimensional Fermi gas is dominated by the Kosterlitz-Thouless tran- 
sition—. 

Throughout this paper, we set fi = k% = 1, and the system volume V is taken 
to be unity. 



2 Two-dimensional trapped Fermi gas in the BCS-BEC crossover 

We consider a two-dimensional Fermi gas, described by the BCS Hamiltonian, 

H = L ^P C po c P° - U L c] p+q/2f-p+q/2l C -p'+^2i- C p'+g/2r C 1 ) 

P a p,p'a 
Here, is a creation operator of a Fermi atom with pseudo-spin a =T,4-, de- 

2 

scribing two atomic hyperfine states. % p — £ p — jx — — /I is the kinetic en- 
ergy measured from the Fermi chemical potential ju, where m is an atomic mass. 
The pairing interaction — U(< 0) is assumed to be tunable. As usual, we mea- 
sure the interaction strength in terms of the s-wave scattering length a s . In the 
two-dimensional case, the scattering length is related to — U a s 13 ' 14 

, m 2 1 

~ u -^ ln 7T^=^— -£ F , .* 0c , ( 2 ) 
2n CyJ—2mE + a s p E + + w— 2e p 
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where 8 = +0, and C = e r (where 7= 0.577 is the Euler's constant). In Eq. (0, 
E+ is a small positive energy, which is taken to be zero in the final expressions for 
the r-matrix theory. 

Figure Q] shows the calculated T c and fx(T = T c ). In panel (a), T c continuously 
changes from the weak-coupling BCS result — r^-Tp to that for a molecular Bose 

f o llamas r 

gas 2 (where Tp is the Fermi temperature). In this BCS-BEC crossover, the 
chemical potential fx at T c changes from /x = £p to — c ^ ma ^ , as shown in panel 
(b). Figure [TJ indicates that the character of superfluidity changes from the BCS- 
type to the BEC of bound molecules around — 1 < In & + 1 ■ Thus, al- 
though there is no clear boundary between the weak-coupling BCS regime and 
the strong-coupling BEC regime, in what follows, we call the region In c £ < 

and In > the BCS side and BEC side, respectively. 

Effects of pairing fluctuations are conveniently described the self-energy cor- 
rection Zp{i(O n ) in the single -particle Green's function, 

G p( ic °n) = ror„ \-i v r77\ ' (3) 

where (O n is the fermion Matsubara frequency, and G p (i(O n ) = (i(Q„ — ^p)" 1 is the 
Green's function for a free Fermi gas. In the T-matrix approximation, the self- 
energy E p (i(0„) is diagrammatically given by the sum of ladder-type diagrams. 
The resulting expression is&i^ 

Z p {i(O n ) = T £ r q (iV n )G° q _ p {iV„ - iOn). (4) 
p.v„ 

Here, V„ is the boson Matsubara frequency. The particle-particle scattering matrix 
r q (iv n ) describes pairing fluctuations in the Cooper channel, given by 

where TI q {iv n ) = T^ P . M , G p+q /2( iv r> + i(0 n) G -p+ q /2^ i(0n ") is the lowest order pair- 
propagator. 

We now include effects of a trap potential. For simplicity, we assume an isotropic 
harmonic trap (with the trap frequency Q) as V(r) = mQ 2 r 2 /2. Within LDA, ef- 
fects of this trap can be conveniently incorporated into the theory by simply replac- 
ing the Fermi chemical potential jj. with ju(r) = fi — V(r). Then, the LDA Green's 
function G p (i(O n ), self-energy Z p (i(0„), and the scattering matrix I^(/V„), depend 
on r through the position-dependent LDA chemical potential jll(r). The chemical 
potential jj. is determined by the LDA equation for the number N of Fermi atoms, 
given by 

N = 2n J rdrn(r), (6) 

where n{r) = TY,pa G po(iO)n) is the number density of Fermi atoms in LDA. 

As in the three-dimensional case, we determine T c from the Thouless criterion, 
which states that the superfluid phase transition occurs when the LDA scattering 
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Fig. 1 (color online)Calculated T c (a) and fi(T = T c ) (b) in the BCS-BEC crossover regime 
of a trapped two-dimensional Fermi gas. The scattering length a s is normalized by the Fermi 
momentum k^. T c and fi(T = T c ) are normalized by the Fermi temperature and the Fermi 
energy £p, respectively. TMA' and 'BCS' show the results based on the T-matrix theory with 
LDA and the mean-field theory, respectively. In panel (b), 'BEC shows the result based on the 
theory valid for the strong-coupling BEC limit. 

matrix has a pole at q = V„ = 0. In LDA, this condition is first satisfied in the 
trap center, so that one may solve f^=o(iV„ = 0,r = 0)~ l = 0to determine T Q . We 
actually solve this T c -equation, together with the number equation (O. We briefly 
note that, because of the presence of a trapped potential, although the system itself 
is a two-dimensional one, a finite T c is obtained within the present combined T- 
matrix approximation with LDA—. 



3 Local density of states and photoemission spectrum 

Once Fig[T]is obtained, the LDA local density of states p((0, r) at T c can be calcu- 
lated from the LDA Green's function G p (z'ft),/i — > ju(r)) as 

p{(0,r) =£A p (fl),r), (7) 
P 

where A p ((0,r) = — n~ l \m[G pa {i(O n — > (0 + iS,r)] is the single-particle spectral 
weight. The photoemission spectrum A p (co)f((0) is also related to the spectral 
weight A p (co,r) as^ 

A p (co)f(co) = |f J drA p (oo + V(r),r)f(co+V(r)), (8) 

where f(co) is the Fermi distribution function, and 7?p is the LDA radius of the gas 
cloud at T = for a free Fermi gas with N atoms, fp is a rf-coupling constant. Since 
the current photoemission experiments do not have enough spatial resolution, the 
spatial average is taken in Eq. (HJ. 

Figure [2] shows the local density of states p((B,r) near the trap center (r = 
0.01/?f) when T = T c 15 . In the weak-coupling BCS regime shown in panel (a), 
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Fig. 2 (color online)Calculated local density of states p(to,r) at r = 0.01i?p when T = T c . (a) 
BCS side (In < 0). (b) BEC side (In > 0). 
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Fig. 3 (color online)Spatial variation of the local density of states p((ti,r) at T c . The interaction 
strength In c £ is set to be (a) -2, (b) 0, and (c) +1. /?p is the LDA radius of a free gas of N 
Fermi atoms at T = 0. 



one already sees a dip structure around 0) = 0. Since the superfluid order param- 
eter vanishes at T c , this is just the pseudogap associated with pairing fluctuations. 
Since T c is close to the mean-field value when In = —2.0 (See FigJTJ), this re- 
sult is quite different from the three-dimensional case, where the pseudogap does 
not appear in such a weak-coupling regime. This is simply because of a stronger 
pairing fluctuations existing in the two-dimensional system. 

This pseudogap structure becomes remarkable, as one increases the interaction 
strength. In the BEC side shown in Fig|2|b), the overall structure is very close to 
the BCS superfluid density of states, characterized by a finite energy gap with two 
coherence peaks at the excitation edges. 

Figure [3] shows the r-dependence of LDOS at T c . In the weak-coupling BCS 
regime shown in panel (a), the pseudogap structure seen in the trap center gradu- 
ally disappears in the outer region of the gas cloud. In addition, the peaks around 
the edges of the pseudogap is suppressed when r/R^ > 0.5. While a similar spatial 
dependence of LDOS is still obtained in the crossover region (Fig|3jb)), FigOfc) 
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(a)ln[2/(Ck F a s )]=-2.0 (b) ln[2/(Ck F a s )]=0.0 (c)ln[2/(Ck F a s )]=1 .0 
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Fig. 4 (color online)Intensity of photoemission spectrum A p (a)f{(0) at T c . (a) BCS regime 
(In c ^ = —2). (b) crossover regime (In c £ = 0). (c) BEC regime (In c f = 1). In each 
panel, the dashed line shows the position of the back-bending curve in the spectrum. The solid 
line shows the free-particle dispersion, a) + fi = £ p . 

shows that a large gap structure remains to r = in the BEC regime. In this 
strong-coupling regime, since the Fermi chemical potential jj. is negatively large 
(See FigQIb).), the system is expected to be well described by a gas of two-body 
bound molecules. Thus, the large gap structure seen in FigOfc) is considered to 
originate from a molecular binding energy. 

Figure [4] shows the photoemission spectra in the BCS-BEC crossover at T c . 
In the weak-coupling BCS regime (panel (a)), the peak position in the spectral 
intensity is close to the free particle dispersion ft) + fi = p 2 / (2m). In the crossover 
region (panel (b)), in addition to this free particle dispersion, one obtains another 
branch in the negative energy region. In the strong coupling BEC regime (panel 
(c)), a large gap structure is seen in the photoemission spectrum, reflecting a large 
binding energy of a bound molecule. 

We note that, in Fig|4fb), while the upper free-particle-like branch mainly 
comes from single-particle excitations around the edge of the gas (where the 
pseudogap effect is almost absent), the lower branch is dominated by the pseudo- 
gapped local density of states around the trap center. This lower branch exhibits 
the so-called back-bending behavior (dashed line in Figj4fb)), which is consis- 
tent with the recent experiments on a quasi-two-dimensional 40 K Fermi gas 10 . We 
briefly note that such a back-bending curve can be already slightly seen in the 
weak-coupling BCS regime (dashed line in Fig|4fa)). 



4 Summary 

To summarize, we have discussed pseudogap phenomena originating from strong 
pairing fluctuations in a trapped two-dimensional Fermi gas. Using a combined 
T-matrix theory with LDA, we have examined how the pseudogap appears in the 
local density of states, as well as the photoemission spectrum, in the entire BCS- 
BEC crossover region at T Q . We showed that the photoemission spectrum splits 
into the upper and lower branches in the crossover region. The lower branch, 
which is deeply related to the pseudo-gapped local density of states around the 
trap center, exhibits a back-bending behavior, which is consistent with the re- 



7 



cent photoemission experiment on a quasi-two-dimensional Fermi gas. Since two- 
dimensional Fermi gases have recently attracted much attention, our results would 
be useful for the further development of ultracold Fermi gas physics. 

Acknowledgements We would like to thank S. Watabe, T. Kashimura, D. Inotani, and R. Hanai 
for fruitful discussions. R.W. was supported by the Japan Society for the Promotion of Science. 
Y. O. was supported by Grant-in-Aid for Scientific research from MEXT in Japan (22540412, 
23104723, 23500056). 



References 

1. S. Giorgini, S. Pitaevskii, and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008), 
and references therein. 

2. I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008), and 
references therein. 

3. Y. Ohashi, and A. Griffin, Phys. Rev. A 67, 063612 (2003). 

4. J. T. Stewart, J. P. Gaebler, and D. S. Jin, Nature 454, 744 (2008). 

5. J. P. Gaebler, J. T. Stewart, T. E. Drake, D. S. Jin, A. Perali, P. Pieri, and G. C. 
Strinati, Nature Phys. 6, 569 (2010). 

6. S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 80, 033613 (2009). 

7. R. Watanabe, S. Tsuchiya, and Y. Ohashi, Phys. Rev. A 82, 043630 (2010), 
and Phys. Rev A 85, 039908(E) (2012). 

8. S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 82, 033629 (2010). 

9. S. Tsuchiya, R. Watanabe, and Y. Ohashi, Phys. Rev. A 84, 043647 (201 1). 

10. M. Feld, B. Frohlich, E. Vogt, M. Koschorreck, and M. Kohl, Nature 480, 75 
(2011). 

11. J. M. Kosterlitz and D. J. Thouless, J. Phys. C6, 1181 (1973) 

12. C. J. Pethick and H. Smith, in Bose-Einstein Condensation in Dilute Gases 
(Cambridge University Press, 2002, UK), Chap.2. 

13. V. Pietila, arXiv: 1203.60541 

14. S. A. Morgan, M. D. Lee, and K. Burnett, Phys. Rev. A 65, 022706 (2002). 

15. At T c , the LDA scattering matrix diverges as r q {iv n = 0,r = 0) ~ q~ 2 
(q — > 0), because the Thouless criterion is satisfied at this temperature. In 
the two-dimensional case, this leads to the diverging behavior of the analytic- 
continued LDA self-energy at r = 0, so that one cannot calculate p(ft), r = 0) 
at T c . We note that, although the number density n(r) also diverges in the trap 
center because of the same origin, the total number N of Fermi atoms in Eq. 
(O converges. 



